Clamping of ferroelectric and ferromagnetic domain walls in conical spiral magnets 
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that sets in at 27K. The interaction between the uni- 
form and spiral components of the conical spiral favors 
the spiral plane (and hence the induced P) orthogonal to 
M, but it cannot constrain the direction of rotation of 
spins in the spiral (sometimes called chirality or handed- 
ness), which determines the sign of P. For a given M, 
the spiral states with opposite handedness and opposite 
polarizations are degenerate and both can be stabilized 



by electric field cooling [111. Il2|. 
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We study the structure of domain walls in multiferroic magnets with the conical spiral ordering. 
We formulate a simple spin model which has a conical spiral ground state in absence of magnetic 
anisotropics. We find a transition from the regime where ferromagnetic and ferroelectric domain 
walls are clamped to the regime where they are decoupled and derive a continuum model describing 
rotation of the spiral plane at the domain wall. The importance of these results for the switching 
phenomena observed in CoCr2 04 is discussed. 
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Introduction: The simultaneous breaking of time- 
reversal and inversion symmetries, which gives rise to 
coupling between electric dipoles and magnetic moments 
in multiferroic materials, is observed in a number of frus- 
trated magnets with unusual spin orders and involves in- 
teresting fundamental physics [1]. The practical motiva- 
tion for studying these materials is the search for ways 
to control magnetic patterns by an applied voltage and 
to manipulate electric polarization by magnetic field. 

The magnetically-induced fiip of the spontaneous elec- 
tric polarization was first observed in nickel boracide Q , 
where the coexistence of a ferroelectric state with an anti- 
ferromagnetic ordering exhibiting a linear magnetoelec- 
tric effect, gives rise to the magnetization, M, linearly 
coupled to the spontaneous electric polarization, P [3|. 
More recently, the 90°-rotation of electric polarization in 
an applied magnetic field was discovered in rare earth 
manganites with a spiral magnetic ordering 0]. The cy- 
cloidal spin spiral breaks inversion symmetry and induces 
electric polarization parallel to the spiral plane and nor- 
mal to the spiral wave vector An applied mag- 
netic field results in the 90°-fiop of the spiral plane with a 
concomitant rotation of P. The giant peak in dielectric 
constant of DyMnOs observed at the polarization fiop 
transition [8^] was recently related to the motion of the 
90° multiferroic domain walls in ac electric field . The 
magnetically-induced polarization fiip was also observed 
in HoMn2 05, which in an applied magnetic field toggles 
between two different multiferroic phases with opposite 
orientations of P [10]. 

Similar polarization reversals were recently found in 
the CoCr204 spinel, which shows a conical spiral order- 
ing and is both ferroelectric and ferromagnetic 0, [l2| . 
What is surprising about these polarization fiips is the 
absence of linear coupling between M and P in CoCr204. 
The magnetization and polarization in this material are 
induced by two different magnetic transitions: the mag- 
netization appears below 95K as a result of the collinear 
ferrimagnetic ordering of Co and Cr spins, while the elec- 
tric polarization is induced by the spiral spin ordering 



In Ref. [11] it was suggested that the coupling between 
the magnetization and polarization in CoCr204 occurs 
at domain walls, where M and P change sign simulta- 
neously. The coexisting ferrimagnetic and spiral orders 
imply the existence of at least two types of domain walls: 
the ferromagnetic domain wall separating domains with 
opposite magnetization and the chiral or handed domain 
wall where the direction of spin rotation and the induced 
electric polarization changes sign. The polarization re- 
versals in CoCr204 can be explained by clamping of fer- 
romagnetic and chiral (ferroelectric) domain walls. 

In this Letter we study the structure of domain walls 
in conical spiral magnets and the clamping of ferromag- 
netic (FM) and ferroelectric (FE) domain walls. The 
straightforward numerical simulation of the domain walls 
in spinels is impossible, as the domain wall width, 
is large compared to the period of the spiral /. Since 
maintaining the relations, a <^ I <^ w <^ where a 
is the distance between neighboring spins and L is the 
system size, is crucial for studying the clamping of the 
FM and FE domain walls, we devised a simple model of 
the conical spiral state, which is amenable to numerical 
simulations. To our knowledge, this is the first simple 
and rigorously soluble model where a conical spiral is 
stabilized without magnetic anisotropy. The first model 
where this physics was found, but only variationally, is 
the cubic spinel with nearest-neighbor AB and BB in- 
teractions, a complex structure with 6 coupled conical 
spirals [l3[ , which incidentally, has provided the basis for 
understanding the observed behavior of several chromites 

iH, [H, HE]. 
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The model: The model consists of two coupled chains 
of classical spins [see Fig. [1]: one with the FM coupling 
Ji? < between neighboring spins and another with the 
competing FM nearest-neighbor interaction Ji < and 
the antiferromagnetic next-nearest-neighbor interaction 
J2 > 0, which for > \ favors a simple (flat) spi- 
ral ordering with the wave vector Q = arccos 
Hamiltonian of the model is 



where 



4J2 • 



The 



n=l 

+ JintSn • Sn], (1) 

where and Si are spins in, respectively, the "ferromag- 
netic" and "spiral" chain, the spins being unit vectors. 
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In some region of coupling con- 



stants it is possible to find a such that the lowest eigen- 
value X{q,a) has two degenerate global minima at = 
and q = Q. In this case a solution for the weak constraint 
is a linear combination of the two degenerate eigenvec- 
tors. The conical spiral state Eq.([2]) belongs to the class 
of such solutions and satisfies the strong constraints. 

The phase diagram of the model Eq.([T]) in the J2-Jint 
plane is shown in Fig.O For sufficiently strong interchain 
interaction Sn||Sn and the system becomes equivalent to 
a single chain showing either the FM or the fiat spiral 
ordering. The vertical line separating these two states is 
4J2 = |Ji| + |Jf|, while the border line between the fiat 
and conical spiral states is given by = 2 + . 
This analytically obtained phase diagram was confirmed 
by numerical simulations of the two-chain model. 

The rotation of spins in the conical and fiat spirals 
induces an electric polarization P ex Q x 63, where Q is 
the spiral wave vector [s], 0, 0] • In the conical spiral state 
the polarization vector is orthogonal to the magnetization 
Miles ml. 



Ji 
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FIG. 1: The system of coupled ferromagnetic and frustrated 
chains that shows a conical spiral ground state. 

When a weak interchain coupling Jint between the two 
chains is turned on, the "ferromagnetic" chain will ac- 
quire a small spiral component, while spins in the "spiral" 
chain will rotate in the plane orthogonal to the magneti- 
zation vector in the "ferromagnetic" chain with a small 
canting in the magnetization direction. In fact, such a 
conical spiral state. 



Sn = s± [cos{Qn - 
Sn = S± [cos{Qn - 



- sm{Qn + (/))e2] + syes, 

- sm{Qn + (j))e2 



S\\e3, 
(2) 

where (61,62,63) is an orthogonal basis and 5y + 5^ = 
5? + 5^ = 1, is the ground state of Eq.([T]) in a wide range 
of the model parameters. 

This can be rigorously proven using the generalized 
Luttinger-Tisza method and the idea of "forced degener- 



acy" p^, [15[ . In this method one minimizes the exchange 



energy, which is a quadratic function of spins [see Eq.([T])], 
replacing the constraint on every site, = = 1 
(the "strong constraints"), by the "generalized weak con- 



straint", En=i K 



L{1 + a ), where a is a 



real constant. If the minimal-energy spin configuration 
for some a also satisfies the strong constraints, it is a 
ground state of Eq.([T]). 

In momentum space the energy ([T]) can be written. 



Ji cos q-\- J2 cos 2q 
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FIG. 2: The phase diagram of the two-chain model Eq.([T]) 
for Jf = —0.75 and Jint < 0. All exchange constants are 
measured in units of |Ji| = 1. 

Domain walls: Next we discuss the structure of domain 
walls in the conical spiral state. In the isotropic spin 
model Eq.([Tj) the width of the FM domain wall is of the 
order of the system size, w ^ in which case the normal 
to the spiral plane follows the slowly rotating uniform 
magnetization, so that P rotates together with M. To 
make the width of the domain wall finite, we add to the 
exchange energy the magnetic anisotropy term on the 
sites of the ferromagnetic chain. 



(4) 
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which for A > favors magnetization along the z axis. 

The width of the FM domain is then w ~ -k^ 

To ensure the presence of a FM domain wah in the 
system, we used the anti-periodic boundary condition 
for spin projection on the easy axis in the ferromagnetic 



chain, 5^+^ 



For ah other spin components we 



used open boundary conditions to prevent the enforced 
commensurabiUty of the conical spiral state by finite size 
effects. Despite the simplicity of the model, its numeri- 
cal simulations are computationally demanding, as frus- 
trated spin interactions result in a rather complex struc- 
ture of the domain wall. We find the lowest-energy spin 
configuration for systems with up to 250 sites in each 
chain using the parallel tempering method 1^, 17 1, which 
makes possible sampling of a large region of the phase 
space with multiple local minima. 

By varying the anisotropy strength A and, hence, the 
domain wall width, we observe a transition from the 
regime where FM and FE domain walls are clamped 
to the regime where they are decoupled. Figures [31(a), 
(b) and (c) show the magnetization along the easy axis 
direction in the "ferromagnetic" chain, 5^, and the y- 
projection of the local polarization. 



^n+1/2 ^ (Sn X Sn+l)^, 



(5) 



induced in the "spiral" chain for two values of the 
anisotropy constant A (we assume that the chains are 
parallel to the x axis) . 

For relatively wide FM domain walls the direction of 
the rotation of spins in the spiral changes sign across the 
domain wall [see Fig. [31(a) and (b)], corresponding to the 
clamping of FM and FE domain walls responsible for the 
switching phenomena in CoCr204. When the thickness 
of the FM domain wall becomes smaller, the spins keep 
rotating in the same direction across the wall, so that the 
induced P has the same sign on both sides and the wall 
is purely ferromagnetic [see Fig. [S^c)]. 

Whether the FM domain wall will induce a FE domain 
wall or not depends on the balance between the inter- 
chain exchange energy, E\^^^ and the energy cost of the 
FE domain wall, £^fe- Expressions for the energies sim- 
plify in the limit \ Jf\^ |^int|, when <C 1 and the FM 
domain wall plays a role of the rotating magnetic field of 
magnitude | Jint| applied to the spiral chain. If, further- 
more, the spiral wave vector, given by sin^ ^ = j^^^ ^ -> 
is small, the energy of the FE domain wall, where the 
spiral plane rotates at the rate ^ ^ Q together with 
the applied field, is 



^FE ^ J2Q^ 



where S\\ = jj^'g] , while the interchain energy is 



l-2Si 



1 (- 



wJ 



-^lllJintl^. 
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FIG. 3: (Color online) The magnetization (red crosses) and 
the polarization Pn+1/2 — (Sn x Sn+i)^ (blue circles) in the 
domain wall. The first three panels show the minimal-energy 
state when a FM domain wall is enforced by boundary condi- 
tions for J2 =0.5, Jf = —4, Jint = — 0.1 and three different 
values of the magnetic ansitropy: (a) A = 0.01, (b) A = 0.04 
(FM+FE domain wall) and (c) A = 0.05 (FM domain wall). 
The other two panels show the ground state when a FE do- 
main wall is enforced by an electric field -\-E(—E) applied 
at the left (right) end of the spiral chain for: (d) J2 = 0.5, 
Jf - -4, Jint - -0.1, and A - 0.01 (FE domain wall) and 
(e) J2 = 1, Jf = -1, Jint = -1.5, and A = 0.01 (FE + FM 
domain wall). All exchange constants are measured in units 
of |Ji| = l. 



For > the energy of the spiral decreases when 
its plane rotates and the FE and FM domain walls are 
clamped. This surprising result is explained by the fact 
that the spiral energy in the momentum space has a local 
maximum at g = 0, so that the rotation of the uniform 
component, S\\es^ results in an energy decrease, which 
for > ^ exceeds the reduction due to the rotation of 
the spiral plane. 

For a nearly fiat spiral, S\\ <C 1, the dimensionless 
parameter determining the domain wall structure is A = 

^ = T^\IW-- for A > 1, E^^.+E,„, < anci 



the FM and FE walls are clamped, while for A ^ 1 the 
domain wall is purely ferromagnetic. The shape of the 
domain wall for the nearly fiat spiral can be found from 
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the continuum hmit of Eq. ([T]) : 

2 



E 



Uir 



dv 



const, (8) 



where the angles and O describe the rotation, respec- 
tively, of the magnetization and the spiral plane around 
an axis orthogonal to i, and v is the dimensionless coordi- 
nate along the chain, such that the domain wall thickness 
is ~ A. The second term in Eq.® is the interaction en- 
ergy between the ferromagnetic and spiral subsystems. 
Minimizing E, one can show that for A 1 the spin 
rotation axis follows closely the direction of the magne- 
tization: Q-e^ \-'^e" {D [see Fig.Hb]. For A < 1 the 
orientation of the spiral plane is only slightly perturbed 
in the vicinity of the FM domain wall: max(|0|) cx A^, 
while S\\ ^ cos changes sign across the domain 
wall [see Fig. |3^]. 




<r-+P 



Fig. [S^a)]. This non-reciprocity, originating from a dif- 
ference in energies of the FM and FE domain walls seems 
to be a general phenomenon and may explain why the 
magnetization reversal in multiferroic GdFeOs results in 
an almost complete reversal of P, while the electrically- 
induced changes in M are very small [l^]- 

Spirals in frustrated magnets have a period of 10- 
20A (/ - 13.4A for CoCr204), while the typical width of 
a FM domain wall is an order of magnitude larger. Fur- 
thermore, in cubic CoCr204 the lowest-order magnetic 
anisotropy is of fourth order, which makes the domain 
wall width even larger, resulting in the perfect clamping 
of the FE and FM domain walls [19]. 

In conclusion, we presented (a) the first simple and rig- 
orously soluble isotropic model with conical ordering and 
(b) numerical and analytical studies of the structure of 
domain walls in a simple model of multiferroic magnets 
with conical spiral ordering. We found a transition from 
the regime where ferromagnetic and ferroelectric domain 
walls are clamped, to the regime where they are decou- 
pled and identified the dimensionless parameter that con- 
trols the transition. The clamping of ferroelectric and 
ferromagnetic domain walls explains the apparent con- 
servation of P X M (the "toroidal moment" |ll| ) and the 
magnetically-induced polarization reversals observed in 
the CoCr204. 

This work was supported by the Thrust II program 
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(FOM). TAK acknowledges helpful discussions with S. 
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FIG. 4: The sketch of the purely FM domain wall in the 
conical spiral magnet, where M changes sign while P does 
not (a), the FM+FE domain wall that reverses both M and 
P (b), and the FE domain wall, where the direction of spin 
rotation in the spiral is reversed, while M is unchanged (c). 

Figures[3l(d) and (e) show the minimal-energy spin con- 
figurations for the case when an electric field ^E{—E), 
applied at the left (right) end of the "spiral" chain, sta- 
bilizes a FE domain wall, while the boundary conditions 
for spins are open. Since spins in the "spiral chain" have 
no anisotropy of their own, the width of the FE domain 
wall is determined by the anisotropy in the FM chain 
and the interaction between the chains. Depending on 
the strength of this effective anisotropy, we find two dif- 
ferent types of FE domain wall: the one where the spiral 
plane flips, as in the FE domain walls shown in Fig. [3] 
and Fig.[3Kb), and the one where the direction of spin ro- 
tation reverses, while the orientation of the spiral plane 
remains unchanged [see Fig. [11(c)] . 

The clamping of FM and FE domain walls is non- 
reciprocal: the FE domain wall shown in Fig. [31(d) does 
not reverse M, even though for the same set of pa- 
rameters the FM domain wall induces the FE one [see 
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